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In  nonlinear  formulation  we  examine  the  problem  of  the  dynamic  stability 
of  cylindrical  shells  made  of  composites,  e.g.,  glass^reinforced  plastic,  un¬ 
der  dynamic  axial  compression.  Ve  assume  as  given  the  law  of  mutual  displace¬ 
ment  of  the  ends  of  the  shell.  Using  the  Bubnov- Galerkin  method,  the  equa¬ 
tions  of  the  nonlinear  theory  of  shells  in  partial  derivatives  are  reduced  to 
ordinary  differential  equations  with  aperiodic  coefficients.  The  equations 
obtained  were  integrated  numerically  using  a  BESM-2M  computer  for  various 
loading  rates  and  shell  parameters  relative  to  certain  brands  of  glass-rein- 
forced  plastics.  Similar  problems  for  isotropic  metal  shells  have  been  exam¬ 
ined  previously  [1,  2].  The  case  of  dynamic  application  of  external  pressure 
to  an  orthotropic  shell  was  examined  in  \3-\  .  J  c<r-.  ■*>  1  ^  . 


Let  us  examine  the  behavior  of  a  closed,  round,  cylindrical,  reinforced- 
plastic  shell  subjected  to  dynamic  axial  compression  (Fig.  1).  We  will  solve 
this  problem  in  the  geometrically  nonlinear  formulation.  Let  us  assume  that 
the  shell  is  attached  by  its  ends  with  frames  whose  points  can  obtain  certain 
radial  displacements  while  the  frames  themselves  remain  curved.  In  view  of 
the  fact  that  the  scatter  of  the  experimental  values  of  the  critical  loads  is 
most  influenced  by  the  initial  irregularities  in  the  shape  of  the  shell,  let 
us  study  the  behavior  of  shells  that  have  initial  camber.  Let  us  use  for  the 
shell  made  of  a  composite  a  model  of  orthotropic  design.  We  will  consider 
that  the  principle  directions  of  rigidity  coincide  with  the  generatrix  of  the 
cylinder  and  the  cross-sectional  arc.  The  elastic  properties  of  orthotropic 
shells  are  characterized  by  moduli  and  E2  in  directions  x  and  y,  by  shear 
modulus  G,  and  by  Poisson  coefficients  and  /i2  »  which  correspond  to  trans¬ 
verse  strain  along  lines  y  and  x.  Let  us  mention  the  familiar  relationship 
^1^2  “  ^2^1  • 
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Fig.  1. 


Let  us  select  the  coordinate  system  such  that  its 
origin  coincides  with  a  point  belonging  to  the  median 
surface.  Coordinate  x  will  be  read  along  the  genera¬ 
trix,  y  -  along  the  arc,  and  z  -  along  the  normal  to 
the  median  surface;  let  us  consider  coordinate  t  aa 
positive  along  the  direction  to  the  center  of  curva¬ 
ture.  Let  us  designate  the  displacements  of  some  point  of  the  median  surface 
in  these  directions  as  u,  v,  v.  Let  us  write  the  equation  of  motion  of  an 
element  of  the  shell,  when  the  deflections  are  comparable  with  the  thickness, 
in  the  presence  of  initial  imperfections  (when  q  -  0)  [2]: 
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Here  w  and  w0  are  the  total  and  initial  deflections;  y  is  the  specific  weight 
of  the  shell  material;  flexural  rigidities  in  the  axial  and  circular  direc- 
tions  Dt  2  -  Eij2h3/12(l  -  n^2)’,  the  derived  rigidity  d3  -  D1/i2  +  2D8  ; 
the  torsional  rigidity  DG  -  <?h3/12;  9  -  the  stress  function  -  is  introduced 
by  the  formulas 
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To  determine  the  stress  function  *  let  us  use  the  equation  of  strain  compati 
bility: 


where 
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We  will  not  take  into  account  the  inertial  forces  corresponding  to  dis¬ 
placements  of  u  and  v  in  the  median  surface.  Thus,  in  this  solution  we  re¬ 
frain  from  studying  the  propagation  of  elastic  waves. 


Let  us  use  the  problem  solution  plan  as  in  [1,  3].  The  total  deflection 
is  approximated  as  follows: 


•  'Infiy  +  f  •In'ax  »in*?y) ,  (3) 

where  a  -  nrx/L;  0  -  n/R ;  a  is  the  number  of  half-waves  along  the  generatrix; 
n  is  the  number  of  waves  around  the  perimeter  of  the  shell.  The  expression 
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for  initial  deflection  is  obtained  from  (3)  by  substituting  f0  for  f. 

Let  us  consider  that  the  shape  of  the  initial  wave  formation,  character¬ 
ized  by  the  parameter  j(>,  is  "in  resonance"  with  wave  formation  of  the  shell 
during  strain,  and  that  the  only  preassigned  parameter  is  the  value  of  f0 . 
Such  an  assumption  intensifies  somewhat  the  influence  of  the  initial  inaccura¬ 
cies  . 


The  boundary  conditions  for  deflection  w  will  be: 

w  -  0;  32w/9x2  -  0  when  x  -  0,  L. 

Substituting  function  (3)  into  the  right  side  of  strain  compatibility 
equation  (2),  we  get  the  integral  of  this  equation  in  the  form 

•-  [Ciip-fr  <!+♦*)  +  <:;/,♦]  co.2a*  +  C,(/*-/2)  (1+  V)  a»2(Jy  + 

+  C,(/*-$  **c«4a*  +  C«(/»-#  **c»4  [C,/i  +  C;  </*-/*)  *]  x 

X  sin  a  x  sin  y  +  [C,  (J*  —  $  ip*  +  C'6  h  ip]  cos  2  a  *  cos  2  }  y  +  C7  (/*  —  /§)  X 
X  $* cos 4 a z  cos 2 £ y  +  C»  (/*  —  fy  $*  co* 2a  x  a*  4  ji  y  +  C»  (/*  —  /J)  \p  sin  3  a  x  sin  S y+ 

(4) 

+  CM(/*  — /J)  tslnoxsin3?y  — py*/2. 

In  the  expression  for  •  we  introduce  the  term  py*/2,  which  corresponds 
to  basic  compressive  stress  p. 

Parameters  Cx  ...  C10  depend  on  the  values  of  a,  p  which  characterize 
wave  formation  of  the  shell. 

Let  us  solve  Eq.  (1)  using  the  Bubnov- Galerkin  method.  As  a  result,  we 
get  equations  which  connect  the  deflection  parameters  with  the  time-changing 
load: 

+  +  y«-‘  <*c/a**)  (5) 

+  (C»—  *-0.  (5a) 

In  these  equations  the  following  dimensionless  parameters  and  designations  are 
introduced: 


€•  —  /•/*:  &  — //*  P  —  P  R/V El  1  —  mmR/hL; 

The  values  of  Ax  . . .  As  and  B1"  ...  B7  '  correspond  to  parameters  £,  A,  which 
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characterize  the  form  of  wave  formation  of  a  shell  and  the  mechanical  proper¬ 
ties  of  the  material.  Solving  the  problem  in  the  first  approximation,  in  Eq. 
(5)  we  drop  the  inertia  term  and  use  for  parameter  the  expression  that  re¬ 
sults  from  Eqs.  (5)  and  (5a).  Relative  shortening  of  the  shell  along  the 
generatrix  is  defined  by  the  formula 
; 

i_ 

L 

Substituting  into  (6)  the  expressions  for  w,  v0  and  <?,  and  converting  to  di¬ 
mensionless  parameters,  we  get 

A +  ,(?-$)  S'/S-Ml^-Co)  Pt,/32.  (7) 

Henceforth  we  will  consider  that  the  value  of  e  changes  in  proportion  to  time 
e  -  c't/L,  where  c'  is  the  rate  of  convergence  of  the  ends,  in  m/s.  Let 
us  turn  to  dimensionless  parameters 

i-StRILh-,  Kt~t'tRlLhAx**ilAx.  (8) 

Dividing  Exp.  (7)  by  Ax  and  substituting  the  obtained  dependence,  and 
also  substituting  (8)  into  (7),  we  get 


f  t.  rfV  (£-*  1  “•  <‘> 


3lC(P-CS)  S*f* 

32  Ax  A 


-  (C  -  C.)  - 


(?•-$*»  *  C(C-W1>  >«.(?- 


- - * - I’0' 


where 


S1-4?(e'r*(*/W  (i/W  1  V^TtVtl  y- 


Pig.  2. 


Equation  (9)  was  integrated  by  the  Runge-Kutta 
method  using  a  BESM-2M  computer,  with  the  following 
starting  data:  f  -  >  df/dt  -  0  with  t  -  0.  Cal¬ 
culations  were  performed  for  shells  with  initial  cam¬ 
ber  f0  -  0.001  relative  to  various  brands  of  rein¬ 
forced  plastics.  During  the  calculation  we  varied 
the  shell  parameters;  the  value  of  £  was  taken  as  2, 
3.5,  and  5;  the  convergence  rate  of  the  ends  was  tak¬ 
en  constant  and  equal  to  c'  -  3.5  m/s.  All  curves 
were  constructed  so  that  parameter  t'  -  tAx /0. 605  - 
-  c'JSV0.605  Lh  was  taken  as  the  base. 
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Figure  2  gives  functions  L(t')  for  shells 
with  parameters  L/R  -  2.2;  h/r  -  1/250  for  various 
ratios  E\/Ez  in  the  case  $  -  2.  As  is  known,  t’  -  1 
corresponds  to  the  convergence  of  the  ends  eu  at 
which  the  upper  critical  pressure  is  reached.  From 
the  figure  it  is  evident  that  the  rapid  increase  in 
deflections  begins  earliest  at  curve  1,  for  which 
E\  ”  &2  “  1.  n  -  14,  then  curve  2  with  -  Ez  “2, 
n  -  15,  and  curve  3  with  Ez  -  Ez  -  5,  n  -  14.  Fig¬ 
ures  3  and  4  give  analogous  curves  for  shells  with 
parameters  L/R  -  3.5  and  2.2;  h/R  -  1/250  when  £x  - 

-  Ea  -  1.  In  Fig.  3  (when  £  -  1)  in  the  case  L/R  - 
-3.5  the  curve  for  which  £  —1.3  deviates  earliest 
from  the  abscissa;  for  the  case  L/R  -  2.2  the  curve 
to  the  left  corresponds  to  n  -  14.  In  Fig.  4  (when 
n  -  1)  the  curve  to  the  left  corresponds  to  n  -  7  in 
the  case  L/R  -  3.5  and  when  n  -  8  in  the  case  L/R  - 

-  2.2.  Figure  5  shows  the  dynamic  coefficient  vs. 
the  ratio  Ex  -  Ez  for  shells  with  parameters  L/R  - 

-  2.2;  h/R  *•  1/250,  n  -  2.  As  can  be  seen  from  the 
figure,  the  lift  of  a  shell  increases  with  increasing 
ratio  Et  -  Ez .  With  an  increase  in  the  parameter  n 
the  number  of  depressions  along  the  peripheral  arc 
decreases.  With  increasing  ratio  R/h  the  dynamic  co¬ 
efficient  increases.  Thus,  the  dynamic  effect  during 
the  axial  loading  of  a  cylindrical  shell  made  of  a 
composite  is  expressed  in  the  appearance  of  higher 
forms  of  stability  loss  and  an  increase  in  the  criti¬ 
cal  load  parameter.  With  an  Increase  in  the  degree 
•f  anisotropy,  corresponding  to  an  increase  in  shell 
rigidity  in  the  longitudinal  direction,  the  dynamic 
coefficient  increases. 
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